Dynamic response of one-dimensional interacting fermions 
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We evaluate the dynamic structure factor S{q, ui) of interacting one-dimensional spinless fermions 
with a nonlinear dispersion relation. The combined effect of the nonlinear dispersion and of the 
interactions leads to new universal features of S{q,u}). The sharp peak S{q,u!) oc q5{Lj — uq), 
characteristic for the Tomonaga-Luttinger model, broadens up; S{q,ui) for a fixed q becomes finite 
at arbitrarily large to. The main spectral weight, however, is confined to a narrow frequency interval 
of the width 5u) ~ q^/m. At the boundaries of this interval the structure factor exhibits power-law 
singularities with exponents depending on the interaction strength and on the wave number q. 
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Low-energy properties of fermionic systems are sensi- 
tive to interactions between fermions. The effect of in- 
teractions is the strongest in one dimension (ID), where 
single-particle correlation functions exhibit power-law 
singularities, in a striking departure from the behavior in 
higher dimensions. Much of our current understanding 
of ID fermions is based on the Tomonaga-Luttinger (TL) 
model 0]. The crucial ingredient of the model is the as- 
sumption of a strictly linear fermionic dispersion relation. 
The TL model, often used in conjunction with a powerful 
bosonization technique , allows one to evaluate various 
correlation functions, such as momentum-resolved 0, 
and local |5| single-particle densities of states. 

Unlike the single-particle correlation functions, the 
two-particle correlation functions of the TL model ex- 
hibit behavior rather compatible with that expected for 
a Fermi liquid with the linear spectrum of quasiparticles. 
For example, the dynamic structure factor (the density- 
density correlation function) 

S{q,Lu) =jdxdte''^'^'-^''\p{x,t)p{0,0)) (1) 

at small q takes the form S'tl(9, w) oc q5{uj — uq). It 
means that the quanta of density waves propagating with 
plasma velocity u are true eigenstates of the TL model; 
these bosonic excitations have an infinite lifetime. 

Below we show that such a simple behavior is an arte- 
fact of the linear spectrum approximation. In reality, the 
spectrum of fermions always has some nonlinearity, 

^n/L.,k=±vk + ky2m + ..., (2) 

where the upper/lower sign corresponds to the right/left 
movers (R/L), and k = P^pf are momenta measured 
from the Fermi points ±pf- (For Galilean-invariant sys- 
tems the expansion (jJl terminates at fc^.) 

The finite curvature (1/m ^ 0) affects drastically the 
functional form of S{q,uj). In a clear deviation from the 
results of TL model, power-law singularities now arise 
not only in the single-particle correlation functions, but 



in the structure factor as well. We show that the singu- 
larities in these two very different objects have a common 
origin, proliferation of low-energy particle- hole pairs, and 
evaluate the corresponding exponents. 

Because of the success of the bosonization tech- 
nique 0, it is tempting to treat the spectrum nonlin- 
earity as a weak interaction between the TL bosons. In- 
deed, the nonlinearity gives rise to a three-boson inter- 
action with the coupling constant cx 1/m 6]. However, 
attempts to treat this interaction perturbatively fail as 
the finite-order contributions to the boson's self-energy 
diverge at the mass shell A reliable method of re- 
summing the corresponding series is yet to be developed. 




FIG. 1: (a) Without interactions, the structure factor differs 
from zero only for lj-<uj < uj+. In the presence of interactions, 
S at uj > U1+ as well. For repulsive interactions, S{q,uj) 
has a power- law divergence at lj = ^-{q) with an exponent 
depending on q. (b) Sketch of S{q, w) at a fixed small q <^ pf, 
see Eqs. ^ and 

We found it more productive to approach the prob- 
lem from the fermionic perspective. One has then a 
benchmark reference point: the structure factor of free 
fermions. At zero temperature, the structure factor has 
a simple physical meaning of the absorption rate of a 
photon with energy u and momentum q Q. Without 
interaction, absorption of a photon results in a creation 
of a single particle-hole pair. For a fixed q < 2pF, the 
energy of this pair lies within the interval 

< uj < uj^, uj± = uq ± q^/2m (3) 

(for free fermions u ^ v). The bounds lo+ corre- 
spond to particlc-holc pairs in which a hole (a particle) 
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is created just below (just above) the Fermi energy. If 
Lo is outside the interval (PJ, the energy and momentum 
conservation laws can not be satisfied and the structure 
factor vanishes. Within this interval, S{q, lo) is indepen- 
dent of w. 



S{q,uj) ^ Soiq) ^ m/q, w_ < w < 



(4) 



see Fig. ^ Accordingly, S{q, ui) at a fixed q exhibits a 
peak which has a "rectangular" shape with the width 
Slu = — — q^/m |2j ■ 

If even a weak interaction is now turned on, its effect 
on the structure factor is dramatic, see Fig.^ The struc- 
ture factor still vanishes below the (renormalized) lower 
absorption "edge" lo — oj^, but above the edge S{q,LLj) 
develops a power-law singularity, 



Siq,i^) 
Soiq) 



Slo 



< u ~ w_ ^ Slo. 



The exponent fi here is a smooth function of q, 

■Kq ' 



(5) 



(6) 



{Vp is the Fourier component of the intra-branch inter- 
action potential, see Eq. H1U|I below). The divergence in 
Eq. lO has the same origin as the familiar X-ray edge 
singularity in metals [lfij | . 

In the presence of interactions, absorption of a photon 
is accompanied by the creation of multiple particle-hole 
pairs, which allows the conservation laws to be satisfied 
at arbitrarily high lo so that S{q, lo) ^ Q ior lo > loj^. 
However, at lo = lo^ the structure factor still exhibits a 
power-law singularity. At \lo — w+j <C Su) we find 



S{q,uj) 
Soiq) 



with 
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5uo 



LO < IjJ+ 



, LO > L0+ 



(7) 



\2ttu ) 



(8) 



Here Uq is the interaction between the right- and left- 
movers, see Eq. (|1(JII below. 

Finally, at high frequencies S{q,Lo) can be evalu- 
ated jl2l | in the second order of perturbation theory in 
the interaction between fermions. 



S{q, lo) — 2^ 



uq 



9 9 ? 

LO'^ — u q 



LO — aj+ ^ Slo. (9) 



Note that cx in Eq. jnj, as expected for a mul- 
tipair contribution to the photon absorption rate 
Therefore, for q <ti pf the high-frequency "tail" yields 



a negligible contribution to the /-sum rule: The main 
spectral weight of S{q^ lo) is still confined within the nar- 
row frequency window Q . At the borders of this interval 
S{q, lo) develops power-law singularities, see Eqs. lO and 
{Tj). The resulting rather peculiar shape of the peak in 
S{q, lo) is sketched in Fig.^b). It is certainly very differ- 
ent from a simple Lorentzian assumed in, e.g., Ref. |l3l| . 
If one were to interpret the finite width of the peak as a 
lifetime of the TL bosons, one would conclude that the 
boson's decay is manifestly non-exponential. Instead, it 
is governed by power laws, indicating strong nonlinearity- 
induced correlations between the TL bosons. 

Equations ©-© represent the main result of this Let- 
ter. We now outline their derivation. 

We describe spinless ID fermions by the Hamiltonian 



H = Ec.fc^Lv-^,, (10) 



a.k 



are 



where a ~ R, L, the single-particle energies £_a.k 

given by Eq. 101, Pa,p = T,k'^l,k-p'^a,k^ L is size 
of the system. We consider short-range interaction, so 
that Vb, Uo are finite and Vq — Vp (x (and similarly for 
Up) for small p <^ pp- (Note that Vp — Up for Galilean- 
invariant systems). 




FIG. 2: (a),(b) First-order correction to S{q,Lj), logarithmi- 
cally divergent at tj ^ oj- -I- and lo lo+ — 0. Each arc de- 
notes a Green function for right-movers. (c),(d) Second-order 
contributions that diverge logarithmically at lj lj+ + 0; 
internal loops correspond to left-movers. 

Let us concentrate on the immediate vicinity of the 
lower absorption edge, lo — > lo-. Divergent corrections 
to S{q,Lo) appear already in the first order in the inter- 
action strength, see the diagrams (a) and (b) in Fig. ^ 
Evaluation of these contributions (which differ from zero 
only at cj- < lo < lo^) yields 



SS{q,io) 
Soiq) 



= p In 



Slo 



< LO - Lo^ ^ Slo (11) 



with p given by Eq. © . It is not difficult to pinpoint the 
origin of the logarithmic divergencies in the perturbation 
theory. Absorption of a photon results in the creation of 
a "deep" hole at fc w —q. Particles near the Fermi level 
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may then scatter off the hole with a small (<C q) momen- 
tum transfer. Excitation of multiple low-energy particle- 
hole pairs then leads to the power-law enhancement of 
the absorption rate, similar to the edge singularity in the 
X-ray absorption spectra in metals [lO|. Unlike the con- 
ventional X-ray singularity problem, in our case the deep 
hole is mobile. It is known, however, that in ID the edge 
singularity remains intact even when the dynamics of the 
hole is taken into account [l^ . 

The above perturbation theory analysis and the anal- 
ogy with the X-ray singularity suggest that the most 
divergent terms of the perturbative expansion can be 
summed by replacing the original model Eq. Hl()|l with an 
appropriate effective Hamiltonian (ID] . It should include 
two narrow (fco ^ q) strips of states: one around the 
(right) Fermi point A: = 0, and another around k = — g. 
The former accommodates low-energy particle-hole pairs, 
while the later hosts a "deep" hole. The corresponding 
effective Hamiltonian iJ_ is then obtained by projecting 
Eq. Hl()|l onto the states of right-movers with ]fc] < fco <C g 
and Jfc -I- g] < fco (the r and d-subbands in Fig. OJ, while 
states outside these intervals are regarded as either empty 
or occupied. 



-q 
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FIG. 3: The states of right-movers included in the effective 
Hamiltonian I12|l . 

Furthermore, for kg <ti q the spectrum within the two 
subbands can be linearized. Using 

\k\<ko ^ \k+q\<ko ^ ^ 

we write the pro jected Hamiltonian in the coordinate rep- 
resentation |l6|. 

i,_./d.*.(-.„„4),v+/^..*;(-^-»-,4)*. 

- (K) - Vg)Jdxpd{x)prix). (12) 

Here pr{x) = : il;l{x)ip^{x) : and pd{x) — 'il;^{x)ip\{x) are 
densities of particles and holes in the corresponding sub- 
bands (the colons denote the normal ordering). The ve- 
locities Mo and U-q are given by 



= V + Vo/2tt +p/r 



(13) 



with p = 0, —q, and include corrections due to both the 
interaction and the spectrum nonlinearity (we neglected 
Vp — Vq (X p^ here). Finally, uj^ is the lower absorption 
edge given by Eq. Q with u — uq. 



In terms of the effective Hamiltonian l)12|l . the struc- 
ture factor Eq. takes the form 

S{q,oj)= [dxdte"^\B{x,t)B'' {0,0)). (14) 



The operator ijt — ipl-ipd creates an inter-subband 
particle- hole pair (an exciton). With the effective Hamil- 
tonian H12|l . the correlation function H14|l can be eval- 
uated by known methods [l^ Il5| . Indeed, the total 
number of d-holes = ^dxpdix) commutes with H_. 
Since the entire d-subband lies below the Fermi level, the 
ground state of corresponds to Nd = 0. The opera- 
tor i?^ in Eq. H14() creates one d-hole, which propagates 
until it is annihilated by the operator B. Therefore, as 
far as the evaluation of Eq. H14|l is concerned, H- can be 
simplified even further by replacing tpdix) — s- Vi}}d{x)V , 
where is a projector onto states with Nd = 0, 1. It is 
easy to see that the projected operators satisfy 

tpd{y)pd{x) = 0, pd(x)ipd{y) = S{x - y)ipd{y), 

from which it follows that [pd{x) , pd{y)] = 0. 

We now bosonize T/i^-field in H12|l according to 



^|Jr{x) ^ ^/koe''^'^''\ [(p{x),ip{y)] ^iTrsgii{x~y), (15) 
which yields i7_ = Hq + SH, where 



Hn = 



uo 
47r 



dx {d.j,(pf 



dx ipU 



,dx)i'd' 



SH = - — {Vo - Vq) dx pd{x) dx(p. 
Ztt 



(16) 



The Hamiltonian can be diagonalized by the unitary 
transformation iJ_ = 

W = 01 dx pd{x)d.j,ip, 



1 Vb-K 



27r u^o 



uo 



'2' 



where p is given by Eq. To the linear order in V, the 
transformation yields = Hq. At the same time, the 
transformation modifies the operator B\ 



B\x) ^ e'^B\x), 



-iW 



fco e 



-i{l-ti./2)ip 



^d. (17) 



Since = Hq is quadratic, evaluation of Eq. (|14() is 
straightforward. The structure factor vanishes identically 
at w < w_, while at w > uj_ it is given by Eq. |SJ|, valid 
with logarithmic accuracy [i.e., up to a numerical factor 
in the square brackets in Eq. 10]. 

Similar procedure can be employed to calculate the 
structure factor near the upper edge lo = w+. At 
uj — !■ UJ+ — the first-order in V correction to S{q, lo) 
[see Fig. nia,b)] diverges. 



SS{q,uj) 
Sa{q) 



= p In 



UJ^ — LU 



5u) 



<uj+- UJ <^5uj. (18) 
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The non- vanishing contributions to S{q, w) at w > uj+ 
appear in the second order in U, see diagrams (c) and 
(d) in Fig. 12 These diagrams describe a process in which 
an absorption of a photon results in the final state that 
has two particle-hole pairs on the opposite branches of 
the Fermi surface jl2j. The corresponding contribution 
reads 

7 / =^ln 1 + )—^: ^>^+ 19) 

oo{q) \ uj - UJ+ J uj + UJ+ 

with 1/ given by Eq. ijSJ. Eq. ((T^ reduces to Eq. © at 
Lu — Lu^ ^ duj, and diverges logarithmically at w — > w+. 

As above, the divergent contributions can be summed 
up by replacing the original model Eq. pU|l with an 
appropriate effective Hamiltonian. In this case the d- 
subband lies well above the Fermi level (near k = q), and 
contains at most a single high-energy particle. However, 
unlike at w — > lu^, the interaction with the left-movers 
now has to be explicitly taken into account {T^. The 
counterpart of Eq. (|12|l then reads 

+ J dx\^{Vo ~ Vq)pdPr + UoiPr+ Pd)pi}, (20) 

where the field operators are defined similarly to ipr^d iu 
Eq. H12(l . and pd,Pr, and pi are the normal-ordered par- 
ticle densities. The structure factor is given by Eq. H14|l 
with the exciton creation operator = V'^'/'r- Similar 
to the above, after bosonizing the r and I subbands, the 
effective Hamiltonian H+ can be diagonalized by means 
of a unitary transformation. A straightforward, although 
lengthy calculation then yields Eq. (TJ for the struc- 
ture factor. Note that at w ^ oj+ the edge singularity 
leads to the suppression of the absorption rate. Techni- 
cally, the difference with Eq. (01 arises because, unlike d- 
holes in Eq. (|12|l , d-particles in Eq. H2()|l move faster than 
the particles near the Fermi level: Uq — uq = q/m > 0. 

The above consideration, leading to power-law singu- 
larities in S{q,uj), is limited to zero temperature T — 0. 
However, as long as temperature remains low, T ^ Soj, 
its main effect is to cut off the power-law singulari- 
ties. This amounts to the replacement \uj — uj±\ — > 
max{7rT, \lu - u!±\] in Eqs. © and iQ. 

In conclusion, we found the dynamic structure factor 
S{q^ Lo) of interacting fermions in one dimension, without 
resorting to the TL model. This allowed us to uncover 
certain universal features in the behavior of S{q, uj). The 
structure factor has a threshold, co = w_(q), stemming 
from kinematic constraints, see Fig^ The divergence of 
S at the threshold, Eq. (0, is characterized by exponent 
/i, which is a smooth function of the wave number q. For 
weak interactions, the explicit form of the function p.{q), 
valid at any q < 2pp, is given in Eq. ©. We believe how- 
ever, that the appearance of the power-law divergency 



in S along the entire boundary lu = LO-{q) is a generic 
feature, not limited to weak interactions. Indeed, the 
methods we employed to arrive at Eq. © rely on small- 
ness of p,{q)^ which is achieved at sufficiently small q at 
any strength of interactions V and U . Also, at q ^ 2pp 
the power-law divergence of S is evident from the con- 
ventional Luttinger liquid theory Remarkably, that 
the origin of the threshold singularity in S can be traced 
back to the physics of Mahan exciton in theory of X- 
ray edge singularity and remains universal at any q. 
The same physics dictates the existence of a power-law 
singularity in S{q,uj) at w = uj+{q), see Eq. 0. 

Besides being of a fundamental interest, the knowledge 
of the structure factor with nonlinear dispersion relation 
is crucial for understanding the Coulomb drag, photo- 
voltaic effect, and other phenomena that owe their ex- 
istence to the particle-hole asymmetry. The developed 
theory is also applicable to the inelastic neutron scatter- 
ing off antiferromagnetic spin chains 0, llfll | . 
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